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The connections between the E(5)— models (the original E(5) using an infinite square well, E(5) — 
f} 4 , E(5) — (3 6 and E(5) — /3 s ), based on particular solutions of the geometrical Bohr Hamiltonian 
with 7-unstable potentials, and the interacting boson model (IBM) are explored. For that purpose, 
the general IBM Hamiltonian for the (7(5) — 0(6) transition line is used and a numerical fit to 
the different E(5) — models energies is performed, later on the obtained wavefunctions are used to 
calculate B(E2) transition rates. It is shown that within the IBM one can reproduce very well all 
these E(5)— models. The agreement is the best for E(5) — /3 4 and reduces when passing through 
E(5) — /3 6 , E(5) — /3 8 and E(5), where the worst agreement is obtained (although still very good 
for a restricted set of lowest lying states). The fitted IBM Hamiltonians correspond to energy 
surfaces close to those expected for the critical point. A phenomenon similar to the quasidynamical 
symmetry is observed. 
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I. INTRODUCTION 

Both, the Bohr-Mottelson (BM) collective model [J, [2, |j| and the interacting boson model (IBM) @, H, H, 0] 
have thoroughly been used to study the same kind of nuclear structure problems. Although very different in their 
formulation, both models present clear relationships. In an approximate way, the IBM can be interpreted as the 
second quantization of the BM shape variables |8[ . More detailed connections between both models were studied 
y—{ . during the eighties by several authors 0, |T(J, HH, Eli [H, EH and, more recently, by Rowe and collaborators [l5| • Both 
' models have three particular cases that can be easily solved and for which a clear correspondence can be done. These 
three cases are: i) the BM anharmonic vibrator and the dynamical symmetry (7(5) IBM limit, ii) the BM 7-unstable 
deformed rotor and the dynamical 0(6) IBM limit, and iii) the BM axial rotor and the dynamical symmetry 0(6) 
IBM limit including Q ■ Q -Q interactions [H|,[I||). Note that although it is traditionally accepted the correspondence 
of the dynamical symmetry SU(3) IBM limit to a submodel of the BM, this fact has never been explicitly probed [HI]. 
Each of these cases are assigned to a particular shape using the Hill- Wheeler variables (/?, 7) [17]: spherical, deformed 



o 

00 
O 



with 7-instability, and axially deformed, respectively. For transitional situations the correspondence between the two 
models is difficult, as Rowe said "what is simple in one model will be complicated when expressed in terms of the 
observables of the other". This situation suggests, for the case of transitional Hamiltonians, to look for the connection 
^ ■ between BM and IBM through numerical studies. 

Among the transitional Hamiltonians, a specially interesting case occurs when it describes a critical point in the 
transition from a given shape to another. In general, for such a situation, where the structure of the system can 
change abruptly by applying a small perturbation, both, the BM and the IBM, have to be solved numerically. 
However, recently Iachello has proposed schematic Bohr Hamiltonians that intend to describe different critical points 
and that can be solved exactly in terms of the zeros of Bessel functions. The first of these models is known as E(5) 
[l8l ]. E(5) is designed to describe the critical point at the transition from spherical to deformed 7-unstable shapes. 
The potential to be used in the differential Bohr equation is assumed to be 7— independent and, for the (3 degree of 
freedom an infinite square well is taken. Similar models were proposed later on by Iachello, called X(5) and Y(5) 
[3 120I ] , to describe the critical points between spherical and axially deformed shapes and between axial and triaxial 
deformed shapes, respectively. All these models give rise to spectra and electromagnetic transition rates that are 
parameter free, up to a scale. In spite of their simplicity, some experimental examples were found [2lL |22|. just after 
the appearance of these models. 

In this work, we concentrate on E(5) and related models. It will be published elsewhere the corresponding study for 
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X(5)— models 23]. The formulation of E(5) attracted immediately attention both experimentally and theoretically. 
Soon after the introduction of the E(5) model, the nucleus 134 Ba w as p roposed by Casten and Zamfir [2l| as a 
realization of it. Other experimental examples proposed are: 104 Ru [24| . 102 Pd [25[ and, 108 Pd (26[. Concerning 
theoretical extensions of E(5), first, Arias [27| proposed a generalization of the E2 operator to be used with the E(5) 
model, then Caprio [28| checked that a substitution of the original infinite well in the (3 variable by a finite one, 
which makes the model not exactly solvable anymore, provides similar results. It showed that the E(5) description 
is "robust in nature" , i. e. the main features of the model remain almost unchanged under strong modification of the 
depth of the potential. Arias and collaborators [29l . [3fJ were the first authors who tried to analyze in a quantitative 
way the connection between the U(5) — 0(6) IBM critical point and the E(5) model. In particular, they established, 
looking to few observables that the IBM, at the critical point, gives results close to E(5) for a small (N « 5) number 
of bosons. However, the IBM results for large N nicely reproduce the spectra and electromagnetic transition rates 
of a Bohr Hamiltonian with a f3 A potential (in the following E(5) — (3 A ). Once more, the model is not analytically 
solvable anymore. Levai and Arias (3l| solved the Bohr equation with a sextic potential with a centrifugal barrier 
arriving to almost closed analytical formulae for the energies and wavefunctions. Immediately after, Bonatsos 
and collaborators explored the possibility of getting numerical solutions for the 7-independent Bohr Hamiltonian with 
potentials of the type /3 2 ™, with n > 1 [33|. These sequences of potentials allow to go from the vibrational limit, n = 1, 
to E(5), n — > 00. In particular, in Ref. [33| spectra and transition rates for the potentials /3 4 ,/? 6 and, /3 8 , are given 
explicitly and compared with the original E(5) (infinite square well potential) case. As mentioned above, all these 
models are produced in the BM scheme and a natural question is to ask for the corresponding equivalence in the IBM. 
Is the IBM able for producing the same spectra and transition rates? If yes, does the IBM Hamiltonian correspond 
to a critical point? This work is intended to answer these questions for the E(5) and related models (— f3 4 , — /3 6 and, 
— /3 8 potentials) and analyze the convergence as a function of the boson number. 

For that purpose, a large set of E(S) and related models results for excitation energies and transition rates are 
taken as reference for numerical fits of the general t/(5) — 0(6) IBM transitional Hamiltonian. This procedure will 
allow to establish the IBM Hamiltonian which best fit the different E(b)— models and their relation with the critical 
points. 

The paper is organized as follows: in section |TT] the fitting procedure is described and the obtained results are 
commented. Section IIIII is devoted to study the energy surfaces of the fitted IBM Hamiltonians and to analyze 
these in relation to the critical point. In section UVl the connection between the present results and the concept of 
quasidynamical symmetry is discussed. Finally, in section|V]the summary and conclusions of this work are presented. 



II. THE IBM FIT TO £(5)-MODELS 



A. The model 



The most general, including up to two-body terms, IBM Hamiltonian can be written in multipolar form as, 



H = e d h d + K a P^P 
where fid is the d boson number operator, and 
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The symbol • stands for the scalar product, defined as Tl-Tl — J2m(~ 4 ) TlmTl-m where Tlm corresponds to the 
M component of the operator Tl- The operator jt m — (—l) m ^i^ m (where 7 refers to s and d bosons) is introduced 
to ensure the correct tensorial character under spatial rotations. 

The electromagnetic transitions can also be analyzed in the framework of the IBM. In particular, in this work we 
will focus on the E2 transitions. The most general E2 transition operator including up to one body terms can be 
written as, 
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where e e // is the boson effective charge and x is a structure parameter. 

The £7(5)— models are intended to be of use for 7-unstable nuclei having 0(5) as symmetry algebra. For the 
construction of an IBM 7-unstable transitional Hamiltonian it is sufficient to impose in Eq. ([1]) n 2 = (this implies 
that no Casimir operator from the SU (3) algebra is included) as can be observed if the Hamiltonian (TT|) is rewritten 
in terms of Casimir operators (the definition for the Casimir operators have been taken from (34|): 

H = ^N(N + 4) + (e d + ^ K4 )C 1 [U(5)] + ^K 4 C 2 [U(5)} 

+ K-^-g)G 2 [0(3)] + (|-l K4 )(7 2 [0(5)]-|(7 2 [0(6)]. (8) 

If additionally, we want to construct an IBM transitional Hamiltonian that preserves the 0(5) symmetry, Casimir 
operators for U(5), 0(6) and 0(5) can be included but not the quadratic 0(3) Casimir operator. This condition, 
translated to the multipolar form language used in Eq. ([T]), leads to the constraint K\ — K3/IO — K4/IA = (see 
Eq. ©). In addition, the structure parameter, \> i n the T E2 operator is usually taken as zero in the standard IBM 
calculations for 7- flat Hamiltonians. In our calculations we will impose k 2 = 0, i.e. the 7 flatness. To make more 
simple the later analysis, we will restrict ourself to the case K4 = 0, leaving as free parameters kq, k\, K3 (plus e d 
that fixes the energy scale). In practice, we do not impose the constraint k,\ — K3/IO = but, as it will be shown, the 
condition will be fulfilled in every fit. 



B. The fitting procedure 

In this section we describe the procedure for getting the IBM Hamiltonian which best fit the different £7(5)— models. 
The x 2 test is used to perform the fitting. The x 2 function is defined in the standard way, 

2 _ 1 (Xiidata) - XjjlBM)) 2 

Ndata ~ N r ~- ^ fr. 



par 



where Ndata is the number of data, from a specific £7(5)-model, to be fitted, N par is the number of parameters used 
in the IBM fit, Xi(data) is an energy level (or a £>(£72) value) taken from a particular £7(5)— model, Xi(IBM) is the 
corresponding calculated IBM value, and er^ is an arbitrary error assigned to each Xi(data). 

In order to perform the fit, we minimize the % 2 function for the energies, using Ed, Kq, k\ and K3 as free parameters 
and k 2 and K4 fixed to zero. For doing this task we use MINUIT [351 ]. which allows to minimize any multi- variable 
function. 

The labels for the energy levels follow the usual notation introduced for the £7(5) model: £ enumerates the zeros 
of the (3 part of the wave function, and r is the label for the 0(5) algebra, i.e. the 0(5) seniority quantum number, 
which is a good quantum number along all the transition from f7 (5) to 0(6). The selected set of levels included in 
the fit for the different E(5)— models are: 

• For the £ = 1 band, all the states with angular momentum lower than 8 and r < 5. An arbitrary a — 0.001 is 
used for these states except for the 2^ state for which a = 0.0001 is used. This latter value allows to normalize 
all the IBM energies to E(2^) = 1. Note that the energy of the state is fixed arbitrarily to 1 (remind that 
the spectrum is calculated up to a global scale factor). 



For the £ = 2 band, all the states with angular momentum lower than 5 and r < 3. An arbitrary a = 0.01 is 
used for these states. 

For the £ = 3 band, just the states with (L = 0, r = 0) and (L = 2, r = 1) are included. An arbitrary a = 1 is 
used for these states. 



With this selection, the number of energy levels included in the fit, Ndata, is equal to 17. Note that the state 0^ is 
not a real data to be reproduced because we are interested just in excitation energies and therefore the ground state 
is naturally fixed to zero in both, E(S)— models and IBM. In TableQ]the states included in the fit are explicitly given. 

Once the IBM Hamiltonian is fixed for each £7(5)— model by fitting the energy levels, the x 2 function for the B(E2) 
values is constructed without any additional fitting. The only parameter in the £72 operator ([7]), e e //, is a global scale 
and is fixed to give B(E2; 2f — > 0^~) = 100 in all cases (remind that the structure parameter x = in the E2 operator 
for the transitional class going from E/(5) to 0(6) studied here). The transitions calculated are enlisted in Table HT1 
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TABLE I: States included in the energy fit. 
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TABLE II: B(E2) transitions to be calculated. 



C. The results 

We have done fits of the IBM Hamiltonian (fT]) parameters, so as to reproduce as well as possible the energies of 
the states given in TableUand generated by the different £(5)-models: E(5) ~ /3 4 , E(5) - f3 6 , E(5) - /3 8 , and E(5). 
Calculations for the four cited E(5)— models as a function of the boson number, N, have been performed. 

As it was mentioned before, n,i and k 4 are set to zero in Eq. J!}, while Sd, kq, k\ ,and K3 are free parameters in a 
X 2 fit to the energy levels produced by the different E(5)— models. In figured] the value of the x 2 for a best fit to the 
different E(5)— models as a function of N is shown. Different type of lines in figure [1] represent the fit to a different 
E(5)— model as stated in the legend box. It is clearly observed that for any N the agreement between the fitted IBM 
and the E(5) — f} A model is excellent and is getting worse for E(b) — /3 6 , E(5) — (3 s , up to reach E(5) which is the 
worst case. In particular x 2 {E(5) — /3 4 ) w x 2 (-^(5))/1000. It is worth noting that these results change slowly with 
the boson number and in all cases, except for E{5) — j3 A for which the agreement is always excellent, the x 2 value is 
an increasing function of N . 

In figure [2] the variation of the parameters fitted in the Hamiltonian are shown. Note that the best fit parameters 
give rise approximately to the cancellation of the quadratic Casimir operator for 0(3), i.e. k\ w K3/IO. This can be 
quantitatively observed in Table [TTTJ 

To have a clearer idea of the degree of agreement between the fitted IBM results with the data from the E(5)— models, 
numerical comparisons are shown in Table HVl for N = 60. This table includes not only the states used in the fit, 
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TABLE III: Parameters of the IBM Hamiltonians used in table ITVl 
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FIG. 1: (color online). \ 2 f° r the IBM fit to the energy levels of the different _E(5)-models, as a function of N. 



but also an extra set of states not included in it. These allow to control the goodness of the obtained fit since they 
are predicted states which, as we can see, have their counterpart in the 23(5)— models. The agreement for 23(5) — /3 4 , 
23(5) — /3 6 , and 23(5) — /3 8 is really remarkable for all the states. In the case of 23(5), only the £ = 1 band is perfectly 
reproduced while for the bands with £ = 2 and £ = 3 the agreement is poor. 

The IBM calculations presented are done with the usual IBM codes and consequently are restricted, due to numerical 
limitations, to N around 100. However, one should note that for the transitional class studied in this work the 0(5) 
seniority is a good quantum number all along the transition. This allows to diagonalize easily matrices corresponding 
to large number of bosons using the procedure described in Ref. (30j and explore the quality of the fits in the large 
N limit. The results of the x 2 fitting for such calculations are presented in figure [3] as a function of TV. Note that 
the curves presented in this figure do not match exactly with the corresponding ones in Fig. CD in the common N 
range. This is because in the case in which the 0(5) symmetry is imposed the x 2 function is constructed with only 
one state, of those appearing in table HI per seniority. Then, the number of states included in the fit is different, which 
results in slightly different values for the \ 2 fitted function. The main conclusion to be extracted from Fig. [3] is that 
only the model 23(5) — /3 4 is exactly (at least for the states considered in this work) reproduced by IBM Hamiltonians 
with 0(5) symmetry in the large N limit. For the rest of models the discrepancy in the IBM fit slowly increases as a 
function of N. 

As a test for the produced wavefunctions with the fitted IBM Hamiltonian, they are used for calculating 232 
transition probabilities, 23(232). The effective charge (scale parameter) in the 232 operator 0, is fixed so as to give 
B(E2; 2^ — ► 0^ ) = 100, thus no free parameters are left in this calculation. For the 23(232) 's calculated (not a fit) a 
X 2 value has been obtained for each 23(5)— model with an arbitrary a = 10. In figure [4] the corresponding x 2 value 
is plotted as a function of N for all the 23(5)— models considered. Figure 2] shows a clear dependence of \ 2 on TV. 
The x 2 value decreases monotonically as N increases for all the 23(5)— models, except for 23(5). In this last case, x 2 
start increasing for N w 20. For N < 20, 23(5) provides the best agreement while 23(5) — /3 4 is the worst. This fact 
changes when N increases, and for N w 75 already 23(5) — /3 4 , 23(5) — /3 6 , and 23(5) — /3 s provide a similar (excellent) 
agreement while the x 2 value for 23(5) is clearly larger. 

For a quantitative comparison, the 23(232) values for the selected transitions with N = 60 are shown in table fVl In 
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FIG. 2: (color online). Values of the fitted IBM parameters (see text) as a function of N. Different panels correspond to the 
fit to the different E(5) -models: a) E(5), b) E(5) - /3 s , c) E{5) - f3 6 , and d) E{5) - /3 4 . 



this table, it is clear the remarkable agreement between the IBM calculations and E(5)— models. Note the At = ±1 
selection rule. Thus, the wave functions produced by the fit to the energy levels are giving roughly the correct B(E2) , s. 
However, it should be noted that the calculated IBM B(E2) values always increase as a function of N. Therefore, 
looking at the transition rates B(E2 : 2+ -> 0+), B(E2 : 4| -> 2+), B{E2 : 2^ -> 0|), B(E2 : 0+ -> 2+), and 
B(E2 : 0+ -> 2+) in tabled one observes that, already for N = 60, the IBM values are larger than those provided by 
E(5) and E(5) — /3 8 . This is also observed for some transitions in the E(5) — f3 6 model, but for none in the £7(5) — /3 4 
model. Thus, one expects for these models to start giving larger \ 2 values from a given TV value on. The IBM results 
are always lower that the E(5) — f} A ones and both are approaching as N increases. 

In view of the excellent agreement between E(5)— models and the IBM, we can state, that it is impossible to 
discriminate, from a experimental point of view, between a E{5)— model and its IBM counterpart. 



III. THE CRITICAL HAMILTONIAN 



One of the most attractive features of the E(5)— models treated in this work is that they are supposed to describe, 
at different approximation levels, the critical point in the transition from spherical to deformed 7-unstable shapes. 
Since they are connected to a given IBM Hamiltonian, as shown in the preceding section, this should correspond to 
the critical point in the transition from U(5) to 0(6) IBM limits, i.e. this Hamiltonian should produce an energy 

surface with (^j^j = 0. Is this the case for the fitted IBM Hamiltonians obtained in the preceding section? 

Before starting with the discussion it is necessary to establish a measure on how close is a given IBM Hamiltonian to 
the critical point. 

An energy surface can be associated to a given IBM Hamiltonian by using the intrinsic state formalism [l(| [HJ 
which introduces the shape variables {j3, 7) in the IBM. To define the intrinsic state one has to consider that the 
dynamical behavior of the system can be approximately described in terms of independent bosons moving in an 
average field [3^|. The ground state of the system is written as a condensate, |c), of bosons that occupy the lowest- 
energy phonon state, T\: 

\ C ) = -L(rt)», (10) 
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8.967 


8.669 


7.607 


7.222 


7.075 


6.895 


6.266 


6.295 


4+ * 

7,6 


2,3 


8.967 


8.669 


7.607 


7.222 


7.075 


6.895 


6.266 


6.295 


o+ * 

°3,2 


2,3 


8.967 


8.669 


7.607 


7.222 


7.075 


6.895 


6.266 


6.295 


6 + 5* 


2,3 


8.967 


8.669 


7.607 


7.222 


7.075 


6.895 


6.266 


6.295 


4+* 


3,2 


12.854 


10.437 


10.274 


9.509 


9.363 


9.007 


8.015 


8.078 


o+* 
Z 9 


3,2 


12.854 


10.437 


10.274 


9.509 


9.363 


9.007 


8.015 


8.078 



TABLE IV: Comparison of energy levels for fitted IBM Hamiltonians, with N = 60, compared with those provided by the E(5)- 
models (see text). The asterisk marks states not included in the fitting procedure. In the states labeled with two sub-indexes, 
the first one corresponds to E(5), while the second to the rest of models. 



where 



rt 



+ (3 cos 7 g£q H — — (3 sin 7 (d\ + d)_ 2 ) 
V2 



(11) 



(3 and 7 are variational parameters related with the shape variables in the geometrical collective model |12j . The 
expectation value of the Hamiltonian fTJ) in the intrinsic state (ITU|) provides the energy surface of the system, 
E(N,/3,j) = (e|-ff|c). This energy surface in terms of the parameters of the Hamiltonian (JXJ) and the shape vari- 
ables can be readily obtained [37j , 



(c\H\c) = 



N(3 2 



(1 + /? 2 , 
N(N- 1) 



n 9 

DKl K? 

4 



(i+/? 2 r 



4 



18 



2 



7 

5 3 



9 

- K4 
5 



•4k 2 ) + 2V2/3 3 k 2 cos(3 7 ) 



(12) 



The shape of the nucleus is defined through the equilibrium value of the deformation parameters, (3 and 7, which 
are obtained minimizing the ground state energy, (c\H\c). A spherical nucleus has a minimum in the energy surface 
at (3 — 0, while a deformed one presents the minimum at a finite value of (3. The parameter 7 represents the departure 
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N 

FIG. 3: (color online). \ 2 value for the IBM fit to the energy levels of the different J5(5)-models, as a function of N (large N 
limit), for an IBM Hamiltonian with 0(5) symmetry (see text). 



from axial symmetry, i.e. 7 = and j — n/6 stand for an axially deformed nucleus, prolate and oblate respectively, 
while any other value corresponds to a triaxial shape. An additional situation appears when the energy surface is 
independent on 7 but presents a minimum in (3, being the nucleus 7-unstable. It should be noted that for a general 
IBM Hamiltonian including up to two body terms the shape is either axially symmetric or 7-unstable. Moreover, the 
Hamiltonians considered in this work correspond always to the 7-unstable situation. 

With the tools described above one can study phase transitions in the IBM Q. First, the parameters that define 
the Hamiltonian are the control parameters and normally are chosen in such a way that only one of them is a variable, 
while the rest remain constant. The deformation parameters /3 and 7 become the order parameters, although in 
our case the only order parameter is j3. Roughly speaking, a phase transition appears when there exists an abrupt 
change in the shape of the system when changing smoothly the control parameter. The phase transitions can be 
classified according to the Ehrenfest classification [38j. First order phase transitions appear when there exists a 
discontinuity in the first derivative of the energy with respect to the control parameter. This discontinuity appears 
when two degenerate minima exist in the energy surface for two values of the order parameter (3. Second order 
phase transitions appear when the second derivative of the energy with respect to the control parameter displays a 
discontinuity. This happens when the energy surface presents a single minimum for [3 — and the surface satisfies the 

condition \^tM) — 0. In a more modern classification, second order phase transitions belongs to the high order 

V a P J 0=0 

or continuous phase transitions [381 ]. 

To determine whether a given Hamiltonian corresponds to a critical point or not, the flatness or the existence of 
two degenerate minima in the energy surface should be investigated. For the case of one parameter IBM Hamiltonian, 
e.g. Consistent Q (CQF) Hamiltonians 39], it is simple to find an analytical expression for the critical control 
parameter in the Hamiltonian. However, for a general IBM Hamiltonian it is necessary to rewrite the energy surface 
in a special way, as the one presented in Ref . [401 ] . There, the authors manage to write the energy surface of a general 
IBM Hamiltonian in terms of two parameters. The authors make use of some concepts from the Catastrophe Theory 
[4l| to define the two essential parameters, (7*1, ra). In terms of these they find expressions for the locus, in the 
essential parameter space, that gives a critical point at the origin in /3, called bifurcation set, and for the locus that 
gives rise to two degenerate minima, called Maxwell set. For the Hamiltonians considered in section |TT] K2 = and 
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«4 = 0, in these cases T2 = and r\ can be written as, 



a 3 -Mo+£/(JV-l) , , 

ri 2oi+e/(JV-l)-a3' 1 j 



where 



E = 


Sd + 6 ki 




1 


Oi = 


4«o, 




1 


a3 = 








u = 


2 ' 



(14) 

Note that in the large N limit, Ed is proportional to N (see figured]) and therefore (fT3)) can be approached by, 

n ~ — ; . (15) 

This expression agrees with the use of an energy surface derived through a Holstein-Primakoff expansion [42j . 

In this language, a critical Hamiltonian corresponds to r% = 0. In figure [5] the values of r% as a function of N for the 
IBM Hamiltonians obtained from the fit are presented for the different E(5)— models studied. For the E(5) model, 
the fitted IBM Hamiltonian produces r\ = for N ps 7. In the case of the E(5) — (3 8 the value r\ = is obtained for 
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IBM 


E(5) - (3 s 
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IBM 
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B{E2 




-of) 


100 


100 


100 


100 


100 


100 


100 


100 


B{E2 


4f 


-2+) 


167.4 


165.2 


173.3 


170.7 


176.6 


173.9 


183.2 


180.5 


B{E2 


6f 




216.9 


215.4 


231.6 


227.0 


239.8 


233.9 


256.4 


248.7 


B{E2 




-2+) 


167.4 


165.2 


173.3 


170.7 


176.6 


173.9 


183.2 


180.5 


B{E2 




-of) 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


B{E2 


4 


-2f) 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


B[E2 




-2+) 


113.6 


112.8 


121.3 


118.9 


125.6 


122.5 


134.3 


130.3 


B{E2 




-4f) 


103.3 


102.6 


110.3 


108.1 


114.2 


111.4 


122.1 


118.4 


B{E2 


3f 


-2f) 


154.9 


153.8 


165.5 


162.1 


171.3 


167.2 


183.1 


177.6 


B(E2 




-4f) 


62.0 


61.5 


66.2 


64.9 


68.5 


66.8 


73.3 


71.1 


B{E2 


of 


-2+) 


216.9 


215.4 


231.6 


227.0 


239.8 


233.9 


256.4 


248.7 


B(E2 


of 


-2f) 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


B(E2 


n 


-0+) 


75.2 


90.2 


91.2 


95.9 


99.0 


99.6 


112.6 


107.9 


B{E2 


4f 


-2+) 


124.3 


152.3 


156.1 


163.5 


172.0 


170.5 


197.9 


186.5 


B(E2 


2f 


-0+) 


65.7 


89.3 


91.6 


97.9 


103.7 


103.4 


126.6 


115.9 


B(E2 


of 


-2f) 


86.8 


81.6 


107.6 


100.8 


119.0 


112.1 


141.8 


135.4 


B(E2 


of 


-2+) 


123.2 


155.0 


178.5 


182.0 


205.3 


198.5 


257.9 


235.6 



TABLE V: B(E2) values obtained, for N = 60, for fitted IBM Hamiltonians (see text) compared with those provided by the 
different E(5)— models. 

N 25, while for E(5) — (3 6 it is obtained for N « 70. For the E(5) — (3 A model it is known that r^ = is reached 
for very large number of bosons [1^, [3(| • 

It is worth to show (see figure [6]) that for all the fitted IBM Hamiltonians, the resulting energy surfaces are quite 
flat in a large interval of N values and, therefore, it is justified to say that the fitted IBM Hamiltonians are very close 
to the critical area. As a consequence, the E(5)— models will be appropriated to describe phase transition regions 
close to the critical point. 

IV. QUASIDYNAMICAL SYMMETRIES 

The concept of quasidynamical symmetry (QDS) was introduced in Refs. [HI, 0, 0, [4^, HfjJ and has been used 
in the study of phase transitions. This concept is very useful for working with Hamiltonians that present as limits 
two dynamical symmetries (depending on the value of a control parameter). In this situation, the system shows the 
tendency to hold onto a given symmetry until the control parameter reach a critical value, passing the system, at this 
moment, onto the other symmetry. The remarkable feature is that the system can present a set of states that behave 
as belonging to irreducible representations (irreps) of the corresponding symmetry group, although in fact, they do 
not belong to a given irrep but to a mixture of them. 

In mathematical terms, a QDS can be defined through the embedded representations [HI]: "If a subset of states of 
a system are in one-to-one correspondence with the states of an irrep of a group G and if all the properties of the 
subset of states associated with obscrvables in the Lie algebra of G (including their relationships to one another but 
not necessarily their relationships to states outside of the subset) are as they would be if the states actually belonged 
to an irrep of the group G, then the subset of states is said to span an embedded representation of G" . Therefore, 
in the case of a QDS, there exist a set of states that behave as belonging to a unique irrep of G, although that is 
only apparent, because they correspond to a superposition of irreps, but all their observables (up to certain degree 
of accuracy) are identical to the ones of states within a given irrep. In summary, the states can be expressed as a 
coherent superposition of irreps that behave as a single one. Note that to show that a QDS exists, one has to fix a 
subset of states and the degree of accuracy for the comparison with the observables of the dynamical symmetry. 

In our comparison between the IBM and the E(5)— models we observe a phenomenon which resembles the QDS, 
i.e. part of the IBM spectrum behaves as having E{ 5)— symmetry, although, indeed they do not have such a symmetry. 
We should emphasize that this is not a real QDS for two reasons: i) E(5)— cases are not dynamical symmetry limits 
of the IBM and ii) the BM and the IBM have different Hilbert spaces. Indeed, it is not possible to define irreps 
in E(5)— models and therefore embedded representations. We will call this situation quasi-critical point symmetry 
(QCPS) g3. 
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FIG. 5: (color online). Values of n (see text for definition) as a function of N for the fitted IBM Hamiltonians. 



In order to study in detail the QCPS one has to fix the degree of accuracy to be demanded to the observables. In 
our study, for the energies an accuracy of 1% for all the states belonging to a given £ is set while for the B(E2) values 
an accuracy of 10% for all the studied intra-band transitions in a given £ is selected. 

Tables llVl and IVl which correspond to N — 60 are analyzed below, 

• E(5): only the states in the £ = 1 band present E(b) QCPS. 

• £7(5) - (3 8 : only the £ = 1 states present E(5) - /3 s QCPS. 

• E(5) - f3 6 : only the £ = 1 states present E(5) - (3 6 QCPS. 

• E(5) - (3 4 : all the studied states, £ = 1, £ = 2 and £ = 3, present E{5) - /5 4 QCPS. 

These results, regarding the energies, can be extended to larger values of N too (see figure [TJ, i.e. the values of the 
energies remain stable when N increases, while for the B(E2) values the observed differences become larger, specially 
in the E(5) case. 



V. SUMMARY AND CONCLUSIONS 



In this paper, we have studied the connection between the E(5)— models and the IBM on the basis of a numerical 
mapping between models. To establish the mapping we have performed a best fit of the general 17(5)— 0(6) transitional 
IBM Hamiltonian to a selected set of energy levels produced by several E(5)— models. Later on, a check to the 
wavefunctions, obtained with the best fit parameters, has been done by calculating relevant B(E2) transition rates. 
All calculations have been done as a function of the number of bosons. Once the best fit IBM Hamiltonians to 
the different E(5)— models are obtained, their energy surfaces are constructed and analyzed with the help of the 
Catastrophe Theory so as to know how close they are to a critical point. Finally, the concept of quasi-critical point 
symmetry is introduced, as similar to the idea of quasidynamical symmetry. 

We have shown that it is possible, in all cases, to establish a one-to-one mapping between the E(5)— models and 
the IBM with a remarkable agreement for both the energies and the B(E2) transition rates. In general, the goodness 
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0.1 0.2 0.3 0.4 0.5 

P (a.u) 

FIG. 6: (color online). IBM energy surfaces as a function of /3, for selected values of n (see text for definition). 

of the fit to the energies is independent on the number of bosons, but the corresponding B(E2) transition rates are 
indeed sensitive to N. This is so specially in the £(5), for which the \ 2 value reaches a minimum for N small (N m 7) 
and from there on increases notably as a function of N. Globally, the best agreement is obtained for the E(5) — (3 4 
Hamiltonian and the worst for the E(5) case. For the case of very large number of bosons and Hamiltonians with 
0(5) symmetry we have confirmed the results of [2!|[3(|, i.e. the only E(5)— model that can be reproduced exactly by 
the IBM is E(5) — /3 4 , corresponding such a Hamiltonian with the critical point of the model (ri = 0). A consequence 
of this excellent agreement is that it is impossible, from a experimental point of view, to discriminate between a 
i?(5)-model and its corresponding IBM Hamiltonian when only few low-lying states are considered (usually the four 
lowest states in the ground state band, plus 0^" and 2g" in the £ = 2 band). 

We have also proved that all the E(5)— models correspond to IBM Hamiltonians very close to the critical area, 
< 0.05. Therefore, one can say that the E(5)— models are appropriate to describe transitional 7— unstable 
regions close to the critical point. 

We have found that the results presented in this paper are consistent with the existence of something similar to a 
quasidynamical symmetry, we call this phenomenon quasi-critical point symmetry. 

Finally, it should be noted that the use of a more general 1/(5) — 0(6) Hamiltonian, e.g. using K4 as free parameter, 
do not change the main conclusions of this work. 
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